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Using the simulation code SELFAS2, we present predictions of the radio signal emitted by ex- 
tensive air showers (EAS) during their development in the atmosphere. The radio emission in the 
MHz range coming from air showers is the superposition of two mechanisms: the variation of the 
transverse current due to the systematic opposite drift of electrons and positrons in the Earth's 
magnetic field and the variation of the charge excess due to the electrons in excess in the shower 
front. In this paper, we stress particularly the effect of the realistic air refractive index on the radio 
signal predicted by SELFAS2. 



INTRODUCTION 



SELFAS2 CONCEPT 



In 1965, Jelley et al |T] and in 1971, Allan [2 showed 
experimentally that the extensive air showers induced 
by cosmic rays were detectable using the radio detec- 
tion technique in the MHz range. With the resurgence of 
this technique in the last decade with the new generation 
of radio detector arrays like CODALEMA [3], LOPES 
[1], RAuger [S] and recently AERA [5], various modern 
theoretical approaches to model the radio signal during 
the air shower development emerged [THlOj. After more 
than ten years of modern theoretical studies, the different 
simulation codes available on the market appear to be in 
quite well agreement [llj . showing that the air shower 
radio emission is quite well understood. 

The electric field emitted by the air shower during its 
development is due to two main mechanisms. The first 
mechanism comes from the time variation of the trans- 
verse current caused by the systematic opposite drift of 
electrons and positrons in the geomagnetic field (the ge- 
omagnetic component, see [31 [T^l US] for experimental 
evidences). The second mechanism comes from the time 
variation of the net negative charge due to the electrons in 
excess in the shower front (the charge excess component, 
see [T4HI6] for experimental evidences). In most cases, 
the geomagnetic component is dominant, except for air 
showers parallel to the magnetic field orientation, where 
the Lorentz mechanism is weak even null. In this partic- 
ular case, the proportion of the charge excess component 
to the total field emitted by the air shower, becomes im- 
portant even dominant. Recently, it has been shown that 
the air refractive index plays an important role on the sig- 
nal amplitude, particularly for antennas located close to 
the shower axis (less than 300 m) [101 E] • 

After a brief overview of the code SELFAS2 [S], we 
will present the latest results predicted by the simulation. 
We will particularly stress and confirm the influence of a 
realistic air refractive index on the radio signal emitted 
by air showers. 



SELFAS2 is based on a microscopic description of the 
shower using the concept of "age" and "shower univer- 
sality" first proposed in [18] to study the longitudinal 
development of purely electromagnetic showers. The im- 
plementation of relevant distributions for secondary elec- 
trons and positrons of the EAS extracted from [TUH21) . 
permits to avoid the heavy use of EAS generators to 
generate air showers in SELFAS2 and makes the simula- 
tion completely autonomous. The particles generated in 
SELFAS2 by Monte Carlo simulation, are tracked along 
their trajectory to compute their individual electric field 
contribution to the total electric field emitted by the air 
shower. With this approach, the characteristics of the 
evolutive spatial density of charge (emissive area) in the 
shower and the systematic drift of electrons and positrons 
due to the geomagnetic field are then naturally taken into 
account in SELFAS2. 

The total electric field 8jS 8j function of the observer 
time t seen at the position x, resulting from the sum- 
mation of all the particles of the shower, is described in 
SELFAS2 by: 
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The subscript "ret" means that the quantities inside the 
bracket [...]rot have to be evaluated at the retarded time 
t^ct related to t by the relation of propagation: 



t = trot + VcS 
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where Ri is the distance between the particle i and the 
observer. The life time of the particle i is described by 
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with ±e the charge of the particle (positron or electron), 
t\ , the retarded instant of creation of the moving parti- 
cle i by sudden acceleration from vl = to Vi(t,-et) and 
^2 , the retarded instant when the particle i stops by sud- 
den deceleration from Vi{tret) to vl = 0. The total track 
length of the particle i contained between t\ and ^2 is 
subdivided in short tracks along which the velocity is 
considered as constant. 
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the line of sight, the electric field emitted by the particle 
i can be rewritten as: 
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AIR REFRACTIVE INDEX AND CERENKOV 
EFFECT 

In the new version of SELFAS2 presented here, the 
influence of the air refractive index (different from the 
unity) in the calculation of the signal emitted by the 
shower, is now taken into account and is described in 
Eqjljand in Eq[2]by the variables rji and rj^e respectively. 
The value of the air refractive index rji considered at the 
instant of emission, depends on the altitude of the par- 
ticle i. The model adopted in SELFAS2 to describe 77 as 
a function of the altitude h (in km) is described in f22^ 
and is given by: 



(4) 



where 770 — 325 and a — 0.1218 km^^. The variation 
of the air refractive index with altitude implies that the 
field propagation from the particle i to the observer at 
ground is not a straight line, however, it was shown in 
P5] that the deviation is negligible. As a consequence, 
the propagation of the field in SELFAS2 is considered 
as a straight line between the antenna and the observer 
position. In order to take into account the variation of the 
air refractive index during the field propagation from the 
particle i to the observer, the effective air refractive index 
which has to be considered in [2] is obtained as follow: 
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where hi is the particle altitude at the instant of emission 
and z°^^, the observer altitude. 

The air refractive index different from the unity implies 
that the denominator in EqJT] can vanish when the angle 
Oi between the line of sight fii and the direction of the 
particle propagation f5i, approaches the Cerenkov angle 
given by: 
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Starting again from the integral expression of the electric 
field given by Eq.(ll) in |9] and using the fact that in the 
Fraunhofer approximation, the path difference between 
the beginning and the end of a short track is given by 
the distance between these two positions projected onto 
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(although we use here the Lorentz gauge, a similar and 
more detailed demonstration is available in [24 , using a 
coulombian gauge). When Oi approaches the Cerenkov 
angle 6'^°'', the numerator and the denominator in Eqjs] 
both vanish. Multiplying and dividing F{t) by 5t, we rec- 
ognize the first derivative of the Heaviside-step function: 
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giving us a limit to estimate the electric field when the 
observer is located on the Cerenkov angle. 



SELFAS2 RESULTS 

To underline the effect of the air refractive index, we 
present in this section the results predicted by SELFAS2 
for a vertical air shower induced by a lO^'^ eV proton, 
in the Auger site configuration (geomagnetic field and 
ground altitude). In Figjl] and Fig[2] we show the sig- 
nal observed by antennas located at different distances 
from the air shower axis, on the east side of the shower 
core. The results obtained with an air refractive index 
fixed to unity (n = 1, dotted lines) are compared to the 
results obtained with a realistic description of the air re- 
fractive index (n = n{h), solid lines). At short axis dis- 
tance (less, than 300 m), the impact of the air refractive 
index on the signal amplitude is not negligible, particu- 
larly around 100 m. The air refractive index influence 
as a function of the distance to the axis d, is not obvious 



is not monotonic 



as we see in Figj3j the ratio ^„(h}^^^ 
when d decreases. By analogy with a point-like source 
moving faster than the speed of light in a medium, the 
non- monotonic lateral distribution shown in Fig [3] can be 
considered as an equivalent of the well know "Cerenkov 
ring" that we can observe in Fig|4j However for the case 
of a charge distribution (as it is the case for the shower 
front) the "ring" observed here is a mix of various effects: 

• field enhancement when the observation angle ap- 
proach 0p°'' (see previous section), which is summed 
up for particles of various energies and directions 
(due to angular dispersion); 

• time compression of the signal due to the air refrac- 
tive index effect on the field propagation (see [T7] 
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FIG. 1: East-west component of the electric field seen by antenna 
located at different distances from the air shower axis. The dotted 
lines correspond to the case of an air refractive index fixed to unity 
(n = 1), the solid line to the realistic case (ra = n(h)) 
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FIG. 2: Frequency spectra of pulses shown in Fig^(same legend). 

for instance), which increases the high frequencies 
amphtude, particularly when the observer is close 
to the Cerenkov angle, as it is shown in Fig[l]and 
Figi 

The air refractive index effect on the radio signal fre- 
quently called " Cerenkov effect" is not trivial and should 
not be compared directly to the classical " Cerenkov radi- 
ation" . Additionally to the " ring structure" of the radio 
signal profile at ground, the east-west asymmetry of the 
radio signal due to the interference between the trans- 
verse current contribution and the charge excess contri- 
bution (see [TMlfTl] for more details), impHes a complex 
structure of the radio signal observed at ground as we 
can see in Fig|4] This shows that caution must be taken 
with 1-dimensional profile with azimuthal invariance, fre- 
quently used for experimental studies. 

The characterization of the "ring effect" could be of a 
great interest because its diameter is directly linked to 
the geometry of the air shower and to the distance be- 
tween the source and the observer. Assuming that the 
depth of the maximum of the radio emission is different 



FIG. 3: Lateral distribution of the absolute field strength pulses 
amplitude, observed for antenna located at the north of the ground 
shower core. The dotted lines correspond to the case of an air 
refractive index fixe to unity (n = 1), the full line to the realistic 
case (n = n{h)). 




FIG. 4: Ground footprint of the east-west polarized radio signal 
deposited by the air shower (obtained using the maximum of the 
signal amplitude in the full band). 

for a shower initiated by a proton from that initiated by 
a iron, we can expect that the diameter of the ring could 
be an observable of the nature of the primary. However, 
the observation of such effect requires to record the radio 
signal up to few hundred MHz, to characterize it experi- 
mentally as it is shown in FigjSj 

CONCLUSION 

With this new version of SELFAS2, we confirmed, the 
importance of a realistic description of the air refractive 
index to compute the radio signal emitted by air show- 
ers. For antennas close to the air shower axis, and more 
particularly close to the Cerenkov angle, the signal is 
strongly modified, showing very high frequency compo- 
nent, up to few GHz. Although the effect of the air re- 
fractive is stronger for frequencies above 100 MHz, its in- 
fluence in the experimental frequency windows frequently 
used (below 100 MHz) by various experiments (as CO- 
DALEMA, AERA, RAuger or LOPES for instance), is 
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FIG. 5; Same as Figjijbut using the 1-100 MHz band (top), the 
100-200 MHz band (center) and the 200-1000 MHz band (bottom). 
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still important and can enhance the signal by a factor of 
more than 1.5. 
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